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We study the S = \ quantum antiferromagnetic XY model on finite triangu- 

> 
\^ . lar lattices with N sites in both longitudinal and transverse magnetic fields. 

We calculate physical quantities in the ground state using a diagonalization 

for spins N < 27, and those at finite temperatures using a quantum transfer 

ov 

Monte Carlo method for N < 24. In the longitudinal magnetic field, the 
long-range chiral order parameter seems to have a finite, nonzero value at low 
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temperatures suggesting the occurrence of a classical umbrella-type phase. 
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In the transverse magnetic field, the 1/3-plateau of the magnetization curve 



appears even at low temperatures, in contrast with the classical model. The 

magnetic field dependences of the order parameters suggest that the chiral- 

ordered, the ferrimagnetic, and the spin flop phases appear successively as the 

magnetic field is increased. The transition temperatures are estimated from 

the peak position of the specific heat, and the phase diagrams are predicted 

in both longitudinal and transverse magnetic fields. 
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I. INTRODUCTION 

The antiferromagnetic XY model on the triangular lattice has attracted much interest. 
The Hamiltonian is 

H = 2Jj2(SfS* + SfS]) - H a Y,Sr, (1) 

where J (> 0) is the exchange integral, H a is an external magnetic field along the a direction, 
and (i, j) takes all the nearest neighbor pairs. In the classical model, Miyashita and ShibaEJ 
found using the Monte Carlo method that a long-range chiral-ordered phase exists at low 
temperatures, although the long-range order of the spins disappears at finite temperatures. 
They also suggested that the transition belongs to the universality class of the Ising model. 
The model in magnetic fields was also studied. In a longitudinal field, all the spins lean to 
the z direction but the chiral-ordered phase still exists (Fig. [I], the umbrella- type phase). On 
the other hand, in a transverse magnetic field, Lee et al. showed that the chiral-ordered, the 
ferrimagnetic, and the spin flop phases appear successively as the magnetic field increases 
(in Figs. ^|(a)-(c)).a In the ferrimagnetic phase, the magnetization has a value of 1/3 of the 
saturated value, and the magnetization curve exhibits a plateau, which is often called as a 
1/3-plateau. However, the ferrimagnetic region of the field shrinks as the temperature cools 
down, and vanishes at T = 0. 

In the S = 1/2 quantum model, it has been discussed whether the classical ordered 
phases are broken down by the quantum fluctuation or not.tttj The ground state properties 
of the model on finite lattices have been studied using a diagonalization method.uH The 
largest lattice treated up to now is that with N = 36,0 and the results suggest that the 
chiral order parameter remains finite for N — ► oo. The properties at finite temperatures 
were also studied by several authors using a quantum transfer Monte Carlo (QTMC),BO a 
high temperature series expansion,!^ and a super-effective-field theory.cl Recently, the authors 
carefully studied an anisotropic Heisenberg model using the QTMC method and showed 
that, when XY anisotropy is large, the peak-height of the specific heat increases with iVJlHl 



Thus they predicted that the quantum fluctuation does not destroy the classical chiral- 
ordered phase in finite temperatures. In presence of magnetic fields, on the contrary, it was 
speculated that the quantum fluctuation brings a different ground state from the classical 
one. Chubukov and Golosov studied the ground state of the XY model in the transverse 
field using the spin wave theory, and suggested that the 1/3-plateau appears at T = 0.1121 
However, the spin wave theory does not always give a correct result in the frustrated model.t!-a 
Moreover, when S = 1/2, the critical field at which the plateau appears becomes negative 
because of the l/S'-expansion.EJ 

In this paper, we study the S — 1/2 quantum XY model in magnetic fields. We consider 
the magnetic structures both at T = and at T ^ 0. We calculate quantities in the ground 
state of finite lattices with N < 27 by the diagonalization method. On the other hand, 
those in finite temperatures are calculated on the lattice with N < 24 using the QTMC 
method. From the result, we suggest the following things. (1) In the longitudinal field, the 
umbrella-type phase occurs up to H z = 3.0. (2) In the transverse field, the 1/3-plateau 
really appears even in the ground state. The chiral-ordered, the ferrimagnetic, and the spin 
flop phases occur as the magnetic field is increased. The phase diagrams in both longitudinal 
and transverse fields are predicted. 

We consider the models in the longitudinal and the transverse fields in Sec. II and Sec. Ill, 
respectively. Summary is given in Sec. IV. We explain the QTMC method in Appendix A. 

II. LONGITUDINAL MAGNETIC FIELD 

In this section, we study the model in a longitudinal field H z . First, we consider the 
magnetization and the chiral order parameter in the ground state. The total magnetization 
along the external field H z is given by 

M z = {l/N)Y.Sl. (2) 

i 

The z component of the chirality of the upright triangle at R is defined as 
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where the relation of i — > j — > k is taken counterclockwise, as is shown in Fig. 0(a). The 
eigenvalues of x z (R-) are ±1/2 and 0. The long-range chiral order parameter is defined as 

x 2 = NsWs^T)^ (R)r - (4) 

where R runs over all the upright triangles on the lattice. In Figs. ^| and [|, we show M z 
and y/x 1 , respectively. Since their size dependences are almost negligible, these values will 
be similar to those for N — > oo. The magnetization increases almost linearly with H z , and 
no anomaly is observed until the saturated field H^^ = 3.0. The chiral order parameter 
decreases monotonically with increasing H z , and vanishes at ff^ aa . = 3.0. These results 
suggest that the chiral-ordered phase changes into an umbrella-type one, and it survives up 
to H^x = 3.0. If it is true, the chiral order parameter at H z ^ will be approximately 
given as \fy^ = ^/x I h z =o cos ' 2 @ = Vx I h z =o{^- ~ (M z ) 2 ), where 9 is the angle between the 
spin and XF-plane. We calculate \p\*' and show it in Fig. || together with \/3p. The two 
results are in good agreement with each other. Thus we may conclude that in the ground 
state the umbrella-type phase occurs at H z ^ 0, and it survives up to -ff^ax = 3-0. 



Next we consider (M z ) and J (x 2 ) f° r a finite temperature, where (• • •) denotes the 
thermal average. All the results presented below have been obtained by the QTMC method 
(see Appnexdix A). The number of the states M used in the QTMC calculation are as follows: 
M = 50 for iV = 9-18, M = 10 for JV = 21, and M = 2 for iV = 24. For every lattice except 
for iV = 24, the set of M states is divided into five subsets, and quantities of interest are 
calculated in every subset. Error bars presented in figures given below only mean deviations 
of the values obtained in different subsets. In the QTMC calculation, since the operators 
Ho = 2JJ2(SfS* + SfS]) and Hi = -H z £ S z are commutable, exp[-(W + Wi)/T] may be 
decomposed into exp[— 7i /T] exp[— TCi/T]. Therefore, we can calculate the quantities for a 
finite temperature combining the T transfer and the H z transfer. That is, those quantities 



are transferred from a high temperature, T , to a temperature of interest, Tf in , by operating 
exp[— 7i /T ]. Then, the field is increased by the operator exp[— TCi/Tf in }. In Figs. ^| and [7|, 



we show (M z ) and J (x 2 ) at T = 0.2 for different N, respectively. In low fields, the results 
are very similar to those of the ground state and their size dependences are small. As H z 
is increased beyond H z = 2.5, their size dependences become considerable. Therefore, we 
suggest that at T = 0.2 the umbrella-type phase survives at least H z < 2.5. 



We also consider the temperature dependence of the specific heat C, (M 2 ), and J(x 2 } 
at fixed fields. The specific heat is calculated from 

C = ^r 2 ((H 2 ) - (H) 2 ). (5) 

We found considerable differences in those quantities between for odd N and for even N, but 
they are reduced with increasing N. In Fig. |](a), we plot C as a function of T at H z = 1.0. 
As iV increases, the peak of C becomes higher and sharper at T ~ 0.4. A similar size 
dependence can be seen up to H z ~ 2.0. Results at H z = 2.0 are also shown in Fig. §(b), 
in which the increase of C is seen at low temperatures. In Figs. |S] and |TDL we plot (M z ) and 



(x 2 ) as functions of T. As the temperature is lowered, (M z ) increases slowly and saturates 



at T ~ 0.4. J (x 2 ) increases and its size dependence becomes smaller, which indicates that 



(x 2 ) remains finite in the thermodynamic limit. Therefore all the results suggest that the 
phase transition occurs at least H z < 2.0. 

The phase diagram of the model predicted from the above results is shown in Fig. [TT| . 
The transition line is estimated from the peak position of C. The chiral-ordered phase 
with the longitudinal magnetization (the umbrella-type phase) survives in the longitudinal 
magnetic field H z . 

III. TRANSVERSE MAGNETIC FIELD 

Next we study the model in the transverse field H x . Since the total magnetization M z is 
not good quantum number, we must treat all the 2 N states in the diagonalization method. 



The QTMC method prepares all the Ising states, and it does not suffer from this difficulty. 
We calculate the magnetization, the order parameter defined below, and the specific heat. 
The number of the states M are as follows: M = 50 for N < 18, M = 10 for iV = 21, and 
M = 2 for iV = 24. 

We calculate the total magnetization given by 

M x = J2 S*/N. (6) 

This quantity becomes 1/6 if the ferrimagnetic phase appears and 1/2 if all the spins align 
along the x direction. In Figs. |T2|(a) and (b), we show the magnetization as functions of H x 
at T = 0.133 and 0.333, respectively. The data for different N lie almost on the same line, 
and the temperature dependence is negligibly small. Note that we also calculate (M x ) in 
the ground state for N = 12, and find the result is quite similar to that at T = 0.133. There 
is the 1/3-plateau in the region of 2.0 ^ H x < 4.0. We estimate two edges of the plateau 
H c i and H C 2 at various temperatures, and plot them in the phase diagram shown below. 
These edges depend little on the temperature. Note that even at T = 0.133, (M x ) slightly 
increases with H x between H c i and H c2 , in contrast with that of the Heisenberg model,H3li3 
because (M x ) is not good quantum number in the XY model. 

Now we consider the spin structure. In Fig. 0, we show three possible spin structures 
(a), (b), and (c) which appear in the classical model at low, intermediate, and high fields, 
respectively. These structures are called as the chiral-ordered, the ferrimagnetic, and the spin 
flop phases, respectively. In order to see which structure is realized with H x , we calculate 
the chiral order parameter and x and y components of the sublattice order parameters which 
are defined as 

M a = [\({M%-M a B f + {M% - M«f 

+ {M^-M%) 2 )} 1 ' 2 (7) 

where M"(= X^ec &? /(N/3)) means the a component of the magnetization of the £ sublat- 
tice. In Table. | we show the relation between the order parameters and the spin structures. 



In Figs. |T3| - |l5| , we plot J {x 2 )i M x \ and A4 y , respectively, as functions of H x at T = 0.133. 



When H x < H&, their size dependences are not so large. Especially, for H x < 1.0, \/{x 2 ) 



depends very little on H x . As H x reaches H c i, J (x 2 ) and A4 y decrease markedly and their 
size dependences become larger. On the contrary, A4 X increases and its size dependence 



becomes small. When H c \ < H x < H c2 , J (x 2 ) an d M. v remarkably depend on the system 
sizes. On the other hand, the size dependence of hA x is negligible. When H x > H c2 , A4 X 
decreases and its size dependence is considerable. M. y increases again and its size depen- 



dence is smaller, but J (x 2 ) still decreases. We tentatively assume that the data fit on the 
1/yN-i unction, and estimated the values for N — > oo which are also shown in the figures. 



When H x < H c i, \/{x 2 ), M x : , and M? have a finite, nonzero value. When H c \ < H x < H c2 , 



(x 2 ), and A4 y seem to vanish for iV — > oo, whereas A4 X remains finite. When H c2 < H x , 



A4 X and M y have a finite, nonzero values, but J \x 2 ) disappears. From these results, we 
suggest that the structures (a), (b), and (c) in Fig. ^| occur when H x < H c i, H c i < H X < H c2 , 
and H c2 < H x , respectively. Of course, these picture of the spin structure are compatible 
with the magnetization curve which are shown in Figs. |T^| 

In Figs. [l6|(a)-(d), we present the specific heat C at H x = 1.0, 3.0, 5.0, and 7.0. For 
H x < 5.0, the peak height becomes higher with increasing iV suggesting the occurrence 
of the phase transition. At H x = 7.0, the peak height does not increase with N, but C 
does at low temperatures. Therefore we expect that the phase transition also takes place 
at H x = 7.0. The transition temperatures are estimated from the size dependence of the 
peak-temperature T m x3 



In Fig. [Ty|, we show the phase diagram of the model. This phase diagram is analogous to 
that of the classical model. However, it should be emphasized that the 1/3-plateau widely 
survives at the low temperatures. 



IV. SUMMARY 

In this paper, we have studied the S — 1/2 quantum antiferromagnetic XY model 
on finite triangular lattice in both longitudinal and transverse magnetic fields using the 
diagonalization and the QTMC methods. We have calculated the magnetization, the specific 
heat, and the order parameters for different sizes of the lattices, and examined their size 
dependences to see whether some long-range ordered phase occurs or not. Our results are 
summarized as follows. 

In the longitudinal magnetic field, the long-range chiral order parameter seems to have 
a finite, nonzero value at low temperatures. The peak height of the specific heat increases 
with iV suggesting the occurrence of the phase transition. From the results, we suggest that 
the umbrella-type phase survives at low temperatures, and predict the phase diagram. 

In the transverse magnetic field, the 1/3-plateau is seen in the magnetization curve at 
low temperatures, in contrast with the classical model. The behavior of the order param- 
eters suggests that the chiral-ordered, the ferrimagnetic, and the spin flop phases appear 
successively as the field is increased. The peak height of the specific heat increases with TV 
suggesting the occurrence of the long-range order. Thus we suggest that the three phases 
occur in the transverse magnetic field at finite temperatures as in the classical model. The 
quantum fluctuation destroys the order of the y component of the spin, and the ferrimagnetic 
phase appears at low temperatures. 
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APPENDIX A: THE QTMC METHOD 

In QTMC method, tH<t3 we calculate the following quantity: 

M M 

(A) =Y,ii>k\Aexp(-[3H)\i; k )/Y,{i>k\eM-PH)M (Al) 

k k 

where A is some physical operator and the sum runs over M states each of which is given 

by 

[T 2N 

l^> = VM^ afe|i) ' (A2) 

here C^ is a random number of — 1 < C$. < 1. We can readily show that (A) — » (A) for 

M 

M — > cxd, because (6/M) y^ CjkCjk = Sij + 0(l/vM)x3 A great advantage of using the 

k 

formula (^1) is that we can obtain an approximate value of the average (A) by summing 



up only M terms instead of all the 2 N terms of the Ising states. Using this formula, we 
can treat systems much larger than those treatable by the rigorous formula. Of course, 
statistical errors of 0(l/vM) arise, but we may reduce them as M is increased. 
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FIGURES 
FIG. 1. A spin structure of the classical model in the longitudinal magnetic field, which is 

called by the umbrella-type phase. 

FIG. 2. Three spin structures of the classical model in the transverse magnetic field, (a) The 
chiral-ordered phase, (b) the ferrimagnetic phase, and (c) the spin flop phase. 

FIG. 3. The field dependence of the magnetization M z at T = 0. 

FIG. 4. The field dependence of the chiral order parameter yx at T = 0. 

FIG. 5. 73? and yfj? at T = 0. 

FIG. 6. The field dependence of (M z ) at T = 0.2. 

FIG. 7. The field dependence of y/Jx 1 ) at T = 0.2. 

FIG. 8. The temperature dependence of the specific heat C for different H z : (a) H z = 1.0 and 
(b) H z = 2.0. 

FIG. 9. The temperature dependence of (M z ). 

FIG. 10. The temperature dependence of \J (x 2 ) f° r different H z : (a) H z = 1.0 and 
(b) H z = 2.0. 

FIG. 11. Phase diagram in longitudinal fields. The symbol x is estimated from the peak 
position of the specific heat. The line is a guide to the eye. 

FIG. 12. The field dependence of the magnetization (M x ) for different T: (a) T = 0.133 and 
(b) T = 0.333. 
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FIG. 13. The chiral order parameter yOw vs - H x . The symbol • denotes the extrapolated 
value. 

FIG. 14. The order parameter A4 X vs. H x . The symbol • denotes the extrapolated value. 

FIG. 15. The order parameter M y vs. H x . The symbol • denotes the extrapolated value. 

FIG. 16. The temperature dependence of the specific heat C for different H x : (a) H x = 1.0, 
(b) H x = 3.0, (c) H x = 5.0, and (d) H x = 7.0. 

FIG. 17. Phase diagram in transverse fields. The symbols x and A are estimated from the 
edges of the plateau of the magnetization curve and the peak height of the specific heat, respectively. 
The lines are guides to the eye. 



12 



TABLES 



Structure (a) Structure (b) Structure (c) 

VW) /o =o =o 

M x /O /O /O 

M y /O =0 /O 

TABLE I. The relation between the order parameters and the classical spin structures. 
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